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Abstract 

The optical responses of structured array of noble-metal nanoparticle dimers immersed in a glass matrix 
are investigated theoretically, motivated by the recent experimental observation of the splitting of the 
surface plasmon bands in silver arrays. To capture the strong electromagnetic coupling between the two 
approaching particles in a silver dimer, the spectral representation of the multiple image formula has been 
used, and a semiclassical description of the silver dielectric function is adopted from the literature. The 
splitting of plasmon resonance band of the incident longitudinal and transverse polarized light is found 
to be strongly dependent on the particle diameter and their separation. Our results are shown in accord 
with the recent experimental observation. Moreover, a large redshift for the longitudinal polarization can 
be reproduced. The reflectivity spectrum is further calculated for a dilute suspension of dimer arrays. 
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I. INTRODUCTION 



The optical properties of small metal clusters embedded in a dielectric medium have attracted 
extensive attention in recent years . 1 i 2 i 3 i 4 i 5 The studies have been developed into many new appli- 
cations in nanostructure enhanced spectroscopies such as surface-enhanced Raman scattering and 
single-target molecule detection,^* 7 - with near-field observation on length scales smaller than the 
wavelength of light. It is known that strong absorption of light occurs at certain frequencies due to 
the collective motions of the conduction electrons in metal called the surface plasmon resonance, as 
well as to the inter-band transition of electrons in deeper levels. The plasmon resonant frequency 
depends strongly on the size and the geometry of the particles,— i 8 i 9 i 10 as well as on the dielectric 
function of metal and the surrounding medium.— The studies are significant theoretically because 
these factors lead to characteristic charge confinement and strong local field enhancement.— On 
the other hand, these properties are also of practical importance in the context of future electronic 
and optical device applications. 

For isolated spherical particles with a diameter much smaller than the wavelength of light 
(rf <C A), the interactions between light and metal lead to an homogeneous but oscillating po- 
larization in the nanoparticles, resulting in a resonant absorption peak in the visible part of the 
optical spectrum. The plasmon resonances in more complex structures of nanoparticles such as 
concentric spherical particles, which are the spherical analog of planar multilayer heterostructures, 
can be more complicated. These resonances are due to the hybridization of free plasmons, which 
can be pronounced depending on the geometry of the particles . 13114115 For nanoparticle ensembles 
like metal clusters, however, the electromagnetic coupling between neighboring particles shifts the 
plasmon absorption bands . 16117 For instance, a nanoparticle chain can be utilized for building opti- 
cal waveguides in the nanoscale . 18il9i2 ° These structures exhibit lateral mode confinement smaller 
than the optical diffraction limit, which cannot be achieved with conventional waveguides nor with 
other novel technologies such as photonic crystals or plasmonic stripe waveguides. 

In the linear arrays of nanoparticles, the optical response can generally be anisotropic, because 
the interparticle coupling depends on whether incident light is polarized longitudinal or transverse 
to the chain axis. This is also one of the underlying principles of optically dichroic glass. Nowa- 
days, structured nanoparticle array can be easily made by electron-beam lithography 21 or other 
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fabrication techniques.— On the theoretical side, finite-difference time-domain (FDTD) simulation 
methods can accurately reproduce the spectroscopic properties of plasmon waveguides and optical 
pulse propagation in these structures as well.— Nevertheless, theoretical investigations by the full 
solution of Maxwell's equations are complicated due to the coupling between plasmons of different 
modes. Although there are already some fruitful discussions , 23 ! 24 it is intrinsically a many-body 
interaction problem. Generally, two types of electromagnetic interactions between particles can 
be distinguished: near-field coupling and far-field dipolar interaction depending on the range of 
interactions concerned. For particle separation r larger than the wavelength of light, the far- field 
dipolar interactions with a r _1 dependence dominate. Much work has focused on these far-field in- 
teractions between metal nanoparticles and their possible applications in optoelectronic materials. 
However, relatively little is known about the nature and the properties of the near-field interactions 
of closely spaced metal nanoparticles, which is the object of the present work. The present work 
is related to earlier studies of FDTD by Oliva and Gray,-^ experiment and simulation by Su et 
al.^ the finite elements method by Kottmann and Martin,^ as well as the pioneering works of 
the discrete dipole approximation (DDA) method by Hao and Schatza 2 ^, who all discovered that 
the interparticle spacing in the particle dimers is crucial to their properties. 

In this paper, we will use the multiple-image method 2 ^ and Bergman-Milton spectral 
representatio n 30 ! 31 for a dimer of two approaching particles 32 to investigate the optical extinc- 
tion and reflectance spectrum of structured arrays of noble-metal nanoparticles, motivated by the 
recent experimental observation of the splitting of the surface plasmon resonance bands for polar- 
izations in the arrays . 21 ' 22 By taking into account the strong coupling of two approaching particles 
in a dimer, we will show that the redshift as large as 1.5 eV into the near-infrared regime observed 
in experiment can be understood. The resonant peak broadening and splitting for different po- 
larizations can be predicted from our calculations, and the results for the reflectance spectrum are 
also presented. 

The rest of the paper is organized as follows. In the next section, we review the general 
spectral representation of two approaching particles. In Sec. Hill we examine the normal-incidence 
extinction spectrum and reflectance spectrum of a silver dimer array immersed in a glass matrix, 
which is followed by discussion and conclusion in Sec. HYI 
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II. FORMALISM 



First we review some formulae similar to those appeared in Ref. 32, however, in a much different 
context. Considering an isolated spherical metal particle in a lossless isotropic dielectric matrix 
with real permittivity e m , the complex permittivity of the metal particle is e(o>), where us is the 
frequency of the external electric field E , and will be discussed in detail afterwards. In this case, 
the dipole moment induced inside the particle is 

P = le m pd 3 E , (1) 

o 

where 

$ = ' H = £m (2) 

is the dipole factor, which defines the polarizability of the particle against the host and is related 
to extinction coefficient directly, d is the diameter of the particle. To account for the multipolar 
interaction between a pair of particles (i.e., a dimer) with spacing a (center-to-center distance 
r = a + d) , we use the multiple image formula.— When the dimer is subjected to an unpolarized 
field, the average of the total dipole moment of one particle is given by 

p* =p T (cos 2 9) +p L (sm 2 9) = ^(Pt+Pl), (3) 

where 9 is the angle between the dipole moment and the line joining the centers of the two particles. 
Here pL and px are the longitudinal and transverse dipole moment, respectively, 3 - 
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where a satisfies the relation cosh a = (a + d)/d. Now the new dipole factors (Pl and $t) of a 
particle in the dimer can be extracted for the longitudinal and transverse field case, respectively. 
Using the spectral representation,— we have 

J- n 



ft = E: 
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with the complex material parameter 



(6) 



where 



= Ff ) = -ln(n + 1) sinh 3 ae^ 2n+ ^ a , 
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, n I[l-2 e -^)«], s f) = I[l + e -(^]. (7) 

In case of unpolarized field, the averaged dipole factor (3* can be derived directly from Eq. (j3J), 
Eq. (JIJ) and Eq. (jHJ), namely, 

2 „=1 V s - S ™ S-S'n' J 

Eq. (JSJ) (or Eq. (jSJ)) is an exact transformation of the multiple image expression,— and consists of a 
set of discrete poles Sn and Sn , which deviates from 1/3 (pole of an isolated spherical particle). 

(L) (T) 

In particular, the longitudinal and transverse poles Sn and s n shift asymmetrically to different 
sides from 1/3. That is, an unpolarized field can excite all poles at both sides. The poles almost 
collapse to that of an isolated sphere (si and Sn — ► 1/3) if a ^ d, indicating that the multipolar 
interaction is negligible. However, when the two particles approach to each other and finally touch, 
the longitudinal (transverse) poles decrease (increase) far from 1/3. Thus, in this case, one should 
take into account the effect of multipolar interactions [see Fig. 5 in Ref. 32 for details]. 

The complex dielectric function e{u) is crucial to the optical properties of metal-dielectric 
systems.-! For noble metals, it can generally be described by the free electron Drude-Lorentz- 
Sommerfeld model plus an additive complex contribution from interband transition, i.e., e(u) = 
1 + x DS (u) + x /B (co>). A complicated function of the dielectric dispersion of Ag takes the form 

2 N 

e{u) = 1 + Boo - + £ ^— -, (9) 

j uj J 

which could be adopted to approximate the measured dielectric function over a wide wavelength 
range.— In Eq. 0, cij may be negative. The sum over iV Lorentz functions and the constants 
are meant to reproduce the interband and all other non-Drude contributions to the dielectric 
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function. Lorentz functions are chosen because it is known that e(u) must obey the Kramers- 
Kronig relations. However, in the frequency range of interest (1 ~ 4.5 ev),— a modified Drude 
model is easier to describe the dielectric response of Ag: 33 

e{u)=e h — — v - , 10 

uj[uj + ry) 

with plasmon resonant frequency uj p = 1.72 x 10 16 rad/s and with the high-frequency limit dielectric 
function e h = 5.45, static dielectric function e s = 6.18. These values were fitted out to be in good 
correspondence with the literature values . 11 ' 34 And the collision frequency 7 in the material is 
assumed to be related to the particle diameter d around 10 nm by>2^ 

7=* + ^, (ID 

with bulk Fermi velocity Vf = 1.38 x 10 6 m/s, room temperature electron mean free path I = 52 nm. 
For d = 10 nm, 7 = 3.025 x 10 14 ; for d = 5nm, 7 = 5.785 x 10 14 . These results show that the 
mean free path of an electron in a nanoparticle is reduced compared to its bulk value due to 
inelastic collisions with the particle surface. The 7 values are taken in our latter calculations of 
different metal-dielectric systems, and compared to the two experimental samples, within which 
the diameters of Ag nanoparticles were in the span of 5 ~ 15 nm . 22 ' 35 Note that 7 determines the 
linewidth of the resonant peak. In the diameter range under consideration, Eq. (|11|) is safe^ and 
indicates that a smaller particle diameter d leads to a wider resonant peak. We are not intended to 
quantitatively compare with the experimental data of Ref. 22, otherwise, we would be restricted 
to a somewhat more rigid size- dependent dielectric function, for example, as Westcott et alJ& 
Nevertheless, the local dielectric treatment is satisfactory as Hao and Schatz pointed out that the 
significant effects of size- dependent dielectric responses come to appear for particles with diameter 
less than 5 nm.— 

Let us use £\{u) and e<i(u) to denote the real and imaginary part of the dielectric function 
obtained by Eq. (fTUJ) . respectively, that is E(uS) = £i(uj) + i£2(u). Fig. [T] shows £i(uj) and £2(1^) 
versus light wavelength A in the span of 250 ~ 1500 nm (i.e. photon energy around 0.8 ~ 5eV). A 
negative £i(cu) is guaranteed for the proper phase relation between the external field and particle 
polarization. The dielectric function changes slightly when the particle diameter decreases from 
10 nm to 5nm [not shown in Fig. however, the resonant line shape is very sensitive to 7, i.e. the 
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particle diameter [see Sec. Illlj . ci< A ensures that the plasma resonance is in quasi-static regime, 
so phase retardation is negligible, effects of higher multipoles can also be neglected for isolated 
spherical particle, which means that dipole plasmon resonance dominates.— 

III. NUMERICAL RESULTS 

Now we consider an array of silver dimer immersed in a glass matrix of refractive index 1.61, 
with the spacing between the two particles in a dimer being less than their diameter (a ^ d). Any 
two dimers are assumed to be far away enough, so the dimer-dimer interaction can simply be at far- 
field approximation, which is neglected in our calculation for simplicity. In the particle diameter 
regime around 10 nm, dipole absorption contribution dominates the scattering effect, although 
dipole scattering increases and dipole absorption fades away for increasing particle sizes.— So in 
the quasi-static regime, the extinction coefficient of a well-dispersed collection of small particles 
is mainly contributed by absorption, with absorption cross section proportional to c<jlm(/?).i2£ 
The complex value expressions of dipole factor in Eq. (0), Eq. (jSJ) and Eq. (|SJ) lead to different 
resonant peaks at different frequencies. To calculate optical extinction, (3 is taken for well-dispersed 
(isolated) nanoparticles, while Pl (Pt) is adopted for an array of dimers. 

Figure 121 shows the optical extinction spectra of an array of dimers with particle diameter 
d = 5nm, spacing a is 0.5 nm and 1.5 nm, respectively. For comparison, solid curve is plotted for 
the extinction spectrum of isolated silver particles. The surface plasmon resonant peak is located 
around 410 nm (~ 3.0 eV), which is in agreement with the first ion-exchanged sample irradiated by 
1 Mev Xe in Ref. 22. The sample contains randomly dispersed silver nanocrystals and the resonant 
band is polarization independent. Long-dashed curves and medium-dashed curves in Fig. |21 are 
the extinction spectra for the array of dimers with illumination light polarized in the longitudinal 
and transverse direction, respectively. It is clear that the plasmon resonance band for longitudinal 
polarization is redshifted with respect to that of isolated particles, whereas the plasmon resonance 
band for transverse polarization is blueshifted. These were also observed in the experiment.— In 
detail, after the sample was subsequently irradiated by 30 Mev Si with fluences up to 2 x 10 14 /cm 2 , 
clear alignment of Ag nanocrystals was observed along the ion-beam direction. Additionally, farther 
redshift and blueshift occur when decreasing the spacing of the two nanoparticles in a dimer from 
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1.5 nm to 0.5 nm [see Fig. That is, a stronger electromagnetic coupling induces further band 
shifts.^ 

Similarly, an array of dimer with a larger particle diameter d = 10 nm are investigated in Fig. 01 
for different spacing a = 0.5, 1.5, 2.5 and 3.5 nm. Splitting of the resonant peak for both the 
longitudinal and transverse polarized light can also be observed. However, a second resonant band 
appears as the two particles in the dimer approach to each other. The second peak position is 
around 830 nm (close to 1.5 eV) when spacing decreases to a = 0.5 nm [see Fig. 0(a)], this is 
in good agreement with the experimental observation of the second sample in Ref. 22. For this 
sample, growing and more compact alignment of the silver nanocrystals are assumed to happen in 
response to higher Si ion fluences irradiation (up to 1 x 10 15 /cm 2 ), and much higher irradiation 
fluences induces much larger splitting of the resonant band for both the longitudinal and transverse 
polarizations. In fact, all of these are also obtained from Fig. 01 In principle, even in the quasi-static 
regime, there are different causes of multipeak behavior of optical spectra: (1) the splitting of the 
dipole mode owing to nonspherical particle shapes, (2) the excitation of higher multipole modes 
in irregularly shaped clusters as a result of inhomogeneous polarization (In this case, number of 
resonances strongly increases when the section symmetry decreases), (3) the enhanced excitation 
of multipoles due to image interactions for spheres. Multipeak structures can also be produced 
by appropriate £\(uS) spectra, and may be damped away if £2^) is sufficiently large. But as Fig|T] 
shows, the model dielectric function resulting from experiments of silver nanoparticles doesn't 
exhibit these behaviors, so we are confident that the multipeak behavior observed in the sample is 
substantially due to the strongly coupling between the two particles in a dimer {i.e. (3)), because 
no obvious identical irregular shape can be seen for the samples. 22 

The largest shift of resonant peak of the dimer array is shown for the longitudinal and transverse 
field cases in Fig. HJ at d — 10 nm. The figure shows the farthest shifting of the peak positions 
versus spacing a. The results are obtained by first calculating the extinction spectra with different 
a in the range of 0.5 ~ 20 nm, and then finding out the position of resonant peak at the longest 
(shortest) wavelength for the longitudinal (transverse) polarization. Large shifting is obvious only 
when a jd < 1, and this can also be understood from the spectral representation in the insert of 
Fig. 0(b), which demonstrates that the spectral poles collapse to 1/3 when a tends to be larger 
than 10 nm. Within the spacing of a < d, the redshift for the longitudinal polarization (diamond) 
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is obviously stronger than the blueshift for transverse polarization (circle). The large splitting is 
due to both the dipole mode coupling (collective excitation mode) and the excitation of multipole 
modes by image interactions. 

We also note that the main resonant peak is narrower in Fig. E] than that in Fig. EJ which 
is mainly ascribed to the different intrinsic damping efficient concerning the particle diameter. 
The narrowing of the resonant peak for increasing particle diameter was also observed in the 
experiments.— Note that no obvious peak shifting is observed in the experiments for isolated 
particles after their growing. 

It is instructive to investigate the optical extinction properties of the array of dimers for an 
unpolarized light. We take the average of the longitudinal and transverse dipole factor, i.e. Eq. (jSJ) 
to obtain the extinction spectra. Results are shown in Fig. El for the two different cases of particle 
diameter d — 5 nm and 10 nm, respectively. Corresponding poles and residues of n — 1 to 10 in 
the spectral representation are given in the inserts. It can be clearly seen that extinction spectra 
still change when particles aggregate into structured array from a well-dispersed assemble, even 
though the illumination light keeps unpolarized. This is due to the asymmetric coupling in different 
topology of the field distribution inside and in the vicinity of the particles.— Note that the shifting 
tendency of plasmon-resonance band could be related to the pole deviating from 1/3. For specific 
spacing, say, a = 0.5 nm, the pole of n = 1 of the dimer with large particles [see insert in Fig. [31(b)] 
is further away from 1/3 than that with small particles [see insert in Fig. Efa)]. Equivalently, 
there is a corresponding resonant peak appearing [Fig. E^b)] at long wavelength regime for the 
case of large particles (10 nm), whereas no obvious peak appearing [Fig. Efa)] for the case of small 
particles (5nm). Theoretically, the discrete terms in the spectral representation should generate 
a series of resonant bands, but they are generally very close and superpose to each other, so it is 
hard to resolve them. In many cases the effect is broadening of the resonance only, so the plasmon 
modes remain spectrally unseparated. 

For interest, we compare the reflectivity spectra for the different cases. In the dilute limit, the 
effective dielectric function of the composite is given b y 38 ' 39 

e e = e m + 3e m pf3 } (12) 

where p is the volume fraction of the silver particles. Reflectance at normal incidence is R — 
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1(1 — \fe~ e )/{l + \fel)\ 2 .— Fig. Eland Fig. [7| are plotted as the reflectance versus photon energy of 
the illumination light for the different arrays discussed above, with particle diameter d = 5 and 
d = 10 nm respectively. Volume fraction p = 0.1 [see Fig.£Ja), (b) and Fig. Ufa), (b)] andp = 0.01 
[see Fig. Efc), (d) and Fig. [7(c), (d)] are taken in the calculations. We can see slight shift of the 
reflectivity spectra for light with longitudinal and transverse polarization as compared to the case 
of isolated particles [solid curves in Fig. |H1 and Fig. Ej. Reflectance decreasing in low- frequency 
regime is notable when spacing decreases [see Fig. EJa), (c) and Fig. Ufa), (c)]. A large fluctuation 
in response to longitudinally polarized light is observed in the case of the larger particles with 
a small spacing [long-dashed curves in Fig. Ej. Meanwhile, increasing the volume fraction of the 
metal nanoparticles leads to the descending main peak of the reflectivity spectra. But note that 
the order of the absolute value of reflectance is different for the two different volume fraction. 
These are also ascribed to the very strong electromagnetic coupling between the nearby particles 
and the strongly enhanced local field. 



IV. DISCUSSION AND CONCLUSION 



Here a few comments on our results are in order. Firstly, the tail going up of the extinction 
spectra in Ref. 22 is ascribed to the absorption of the matrix itself, see the reference measurement 
for a Ag ion exchanged sample that was irradiated with Si only. It does not show a plasmon 
absorption band and is colorless, because there were no silver nanocrystals formed. We didn't 
take it into account in our calculation of the extinction spectrum. Furthermore, the experimental 
spectra are broader than our results, this is due to particle size distribution and coupling of particles 
in a chain through many-body interaction although we show that two-body interaction contributes 
little when particles are far away with spacing larger than their diameter. 

In the present work, we are concerned with the multipolar interaction between particles, i.e., we 
focus mainly on the dimerization effect. In doing so, we could neglect the non-local effects in the 
present classical treatment. It is believed that the significant effects of size- dependent dielectric 
responses come to appear for particles with diameter less than 5 nm, 28 Although near adjacence 
may enhance nonlocal contribution, our calculation results depend on the spacing-diameter ratio 
cr/d, rather than the absolute value of a. That is, the spacing in use is scaled by the particle size, 
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if regardless of the retardation effect. So we also didn't consider the suppression of the dielectric 
confinement due to quantum penetration effects with very small interparticle gaps. We should 
worry about retardation effects when the particles size becomes large. 

Due to strong coupling in the dimer, the enhanced excitation of multipoles of the electric field 
that occurs in the vicinity of the dimer appears, which is responsible for the excitation of spectrally 
distinct higher order plasmon modes, such as the quadrupole plasmon-peak feature around 900nm. 
Crosstalk between adjacent dimers is expected to have a negligible effect on the optical spectra since 
it takes place via far-field scattering of the individual nanoparticles with a distance dependence of 
r _1 , whereas near-field interactions of adjacent particles in each dimer show a r~ 3 dependence and 
dominate at small distances, which is captured by our multiple image model. 

Maier et al. used a coupled dipole model in a particle chain and predicted a value of 2 : 1 of the 
ratio of peak shifts of the longitudinal collective modes to peak shifts of the transverse collective 
modes, which is smaller than experiment value 2.3 : 1.— We believe that if one takes into account 
the multiple image interactions in the nanoparticle chain, the theoretically calculated value will 
be increased. Work are under progress in this way. 

We can take one step forward to include the nonlinear characteristics of noble metal particles. 
For instance, based on Eq. ([12)1. we could derive the effective third-order nonlinear susceptibility 
and then the nonlinear enhancement may be studied by taking into account multipolar interactions. 
Regarding nonlinear enhancement due to dimerization, we can formulate some equations, based 
on, say, Yuen and Yu.— It is also interesting to apply the present theory to the polydisperse size 
case, in which the two particles have different diameters. 
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Figure Captions 

Figure [J Complex dielectric function of silver particles obtained from Eq. ()10jl . Parameters: 
d = 10 nm, others given in the context. 

Figure [2J Extinction spectra for an array of dimer with particle diameter d = 5 nm, at two 
different spacing a. The polarization of the incident light is longitudinal (long-dashed curve) 
or transverse (medium-dashed curve) to the axis of the dimer. For reference, solid curve is the 
extinction spectra of isolated and well-dispersed particle collection. 

Figure U2 Splitting of the extinction spectra for an array of dimer with particle diameter 
d = 10 nm, at four different spacing a. Others the same as in Fig. |21 

Figure The resonant peak splitting for an dimer array of diameter d = 10 nm as the spacing 
decreases from 20 nm to 0.5 nm. The lines are guides to the eye. 

Figure Q3 Unpolarized light extinction spectrum for an array of dimer with particle diameter 
d = 10 nm and 5nm respectively, at different spacing a. Inserts in each panel are the poles and 
residues of the spectral representation. 

Figure [BJ Normal incidence reflectivity spectra for an array of dimer with particle diameter 
d = 5nm for different spacing a, taking volume fraction p = 0.1 [(a) and (b)] and p = 0.01 [(c) 
and (d)] respectively. 

Figure Normal incidence reflectivity spectra for an array of dimer with particle diameter 
d = 10 nm. Others the same as in Fig. |H1 
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